ANL-HEP-PR-97-84 
CTEQ-710 
MSUHEP-71027 



Photon Pair Production With 
Soft Gluon Resummation 

Q '. In Hadronic Interactions 

On 

o 
Q 



Balazs/ a ) E.L. Berger,^ S. Mrenna/ 6 ' and C.-P. Yuan( a ) 



Argonne, IL 60439, U.S.A. 



(a) Department of Physics and Astronomy, Michigan State University 
East Lansing, MI 48824, U.S.A. 
(b) High Energy Physics Division, Argonne National Laboratory 

ON 

a,: 

Ph' Abstract 

<D . 
r- 1 _ 

The production rate and kinematic distributions of isolated photon pairs pro- 
duced in hadron interactions are studied. The effects of the initial-state multi- 
^ ■ pie soft-gluon emission to the scattering subprocesses qq, qg, and gg — ► are 

resummed with the Collins-Soper-Sterman soft gluon resummation formalism. 
The effects of fragmentation photons from qg — * jq, followed by q —* jX, are 
also studied. The results are compared with data from the Fermilab Tevatron 
collider. A prediction of the production rate and kinematic distributions of the 
diphoton pair in proton-nucleon reactions is also presented. 



1 Introduction 



An increasing amount of prompt diphoton data is becoming available from the 
Tevatron collider and the fixed-target experiments at Fermilab. A comparison of the 
data to the calculation of the diphoton production rate and kinematic distributions 
provides a test of many aspects of perturbative quantum chromodynamics (pQCD). 
Furthermore, understanding the diphoton data is important for new physics searches. 
For example, diphoton production is an irreducible background to the light Higgs 
boson decay mode h — > 77. The next-to-leading order (NLO) cross section for the 
PP ~ * 77A process |]J was shown to describe well the invariant mass distribution of 
the diphoton pair after the leading order (LO) gg — > 77 contribution (from one-loop 
box diagrams) was included 0, |3[]. However, to accurately describe the distribution 
of the transverse momentum of the photon pair and the kinematical correlation of 
the two photons, a calculation has to be performed that includes the effects of initial- 
state multiple soft-gluon emission. In hard scattering processes, the dynamics of the 
multiple soft-gluon radiation is predicted by resummed pQCD |4|j-[|i"3|j. 

In this work, the Collins-Soper-Sterman (CSS) soft gluon resummation formalism, 
developed for Drell-Yan pair (including W and Z boson) production 0, is extended 
to describe the production of photon pairs. This extension is similar to the formalism 
developed for describing the distribution of the leptons from vector boson decays [fLOf] 
because the final state of the diphoton process is also a color singlet state at LO. 
Initial-state multiple soft-gluon emission in the scattering subprocesses qq, qg and 
gg — > 77X is resummed by treating the photon pair 77 similarly to the Drell-Yan 
photon 7*. In addition, there are contributions in which one of the final photons is 
produced through a long-distance fragmentation process. An example is qg — > 75 
followed by the fragmentation of the final state quark q — > jX. An earlier study of 
soft-gluon resummation effects in photon pair production may be found in Ref. (Tj[ . 



The results of this overall calculation are compared with CDF ]TJ| and D0 [|T(J 
data taken at the collider energy y/S = 1.8 TeV. A prediction for the production 
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rate and kinematic distributions of the diphoton pair in proton-nucleon interactions 
at the fixed-target energy v^S" = 31.5 GeV, appropriate for the E706 experiment at 



Fermilab fll7[ , is also presented. 

Section 2 reviews some properties of the fixed order calculation of the production 
rate and kinematics of photon pairs. Section 3 presents the soft gluon resummation 
formalism and its application to diphoton production. The numerical results of this 
study, and a comparison with data, are presented in Sec. 4. Finally, Sec. 5 contains 
a discussion of the results and conclusions. 

2 Diphoton Production at Leading and Next— to- 
Leading Order 

The leading order (LO) subprocesses for diphoton production in hadron interactions 
are of order a 2 em , where a em denotes the electromagnetic coupling strength. There 
are three classes of LO partonic contributions to the reaction hih,2 — > where 
hi and h,2 are hadrons, illustrated in Fig. (la)-(lc). The first (la) is the short- 
distance qq — > 77 subprocess. The second (lb) is the convolution of the short- 
distance qg — > 7g subprocess with the long-distance fragmentation of the final state 
quark q — > jX. This is a LO contribution since the hard scattering is of order a em a s , 
while fragmentation is effectively of order a em /a s . Here, a s denotes the QCD cou- 
pling strength. Class (lb) also includes the subprocess qq —>■ jg convoluted with the 
fragmentation g — > jX. Finally, there are LO contributions (lc) involving subpro- 
cesses like qq — > qq, where both final state quarks fragment q — > jX. The transverse 
momenta of the photons are denoted Pt x and pr 2 , and the transverse momentum of 
the pair is Qt — Vt x + Pt 2 ■ m the absence of transverse momentum carried by the 
incident partons, the LO process (la) provides Qt = 0. With the added assumption 
of collinear final-state fragmentation, (lb) provides Qt = Vt x + Pt 2 — (1 — z )Vt^, 
where photon 2 carries a fraction z of the momentum of the final-state quark. Given 
a lower limit on the magnitude of the transverse momentum p T of each photon, the 
total cross section at LO is finite. 
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The next-to-leading order (NLO) subprocesses for diphoton production are of 
order ot 2 em ot s . One class of one-loop Feynman diagrams (Id) contributes by interfering 
with the tree level diagram (la). Real gluon emission (le) is also present at NLO. The 
subprocess qg — > 77^ contains a singular piece (If) that renormalizes the lower order 
fragmentation (lb) and a piece (lg) that is free of final-state collinear singularities. 
Finally, subprocesses like qq — *> qq'yy contain a regular piece involving photon emission 
convoluted with a fragmentation function (lh) and pieces that renormalize the double 
fragmentation process (lc). The regular 3-body final state contributions from (le), 
(If), and (lg) provide Qt = — Vr v where j represents the final-state quark or gluon. 
The full set of NLO contributions just described is free of final-state singularities, 

and the total integrated cross section at NLO is finite for a finite lower limit on each 
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Higher order calculations in a s improve the accuracy of predictions for total cross 
sections involving quarks or gluons when only one hard scale Q is relevant. For 
h\h 2 — > 77A, this scale can be chosen proportional to the invariant mass of the 
photon pair, Q = M 77 , which is about equal to 2p T for two well separated photons 
in the central rapidity region. For kinematic distributions that depend on more than 
one scale, a NLO calculation may be less reliable. One example is the distribution of 
the transverse momentum of the photon pair, Qt = \Qt\- At fixed Q, the behavior 
for small Qt of the NLO contribution to the differential cross section has the form 



where ao and a± are dimensionless constants, and <tq(Q) is the total LO cross sec- 
tion calculated from the subprocess (la). The structure of Eq. ([j]) indicates that the 
fixed order QCD prediction is reliable when Qt — Q, but becomes less reliable when 
Qt "C Q, where \ia(Q 2 /Q T ) becomes large. In the region Qt <C Q, the photon pair 
is accompanied by soft and/or collinear gluon radiation. To calculate distributions 

like ^ <7 ) reliably in the region Qt <C Q, effects of multiple soft gluon emission must 
dQ T 

be taken into account explicitly P|-||13||. The contributions (le) and (lg) exhibit 
singular behavior that can be tamed by resummation of the effects of initial-state 




(1) 



4 



multiple soft-gluon radiation to all orders in a s . Other contributions that do not 
become singular as Qt — > do not need to be treated in this manner. Fragmentation 
contributions like (lb) are found to be small in magnitude after isolation restrictions 
are imposed on the energy of the hadronic remnant from the fragmentation. There- 
fore, contributions like (lc) and (lh) are ignored in this work. Gluon fragmentation 
to a photon can be ignored, since its magnitude is small. 

The subprocess gg — > 77, involving a quark box diagram, is of order a 2 m a 2 . 
While formally of even higher order than the NLO contributions considered so far, 
this LO gg contribution is enhanced by the size of the gluon parton distribution 
function. Consideration of the order ot 2 m a1 correction leads to resummation of the 
gg subprocess in a manner analogous to the qq resummation. 

3 The Soft Gluon Resummation Formalism 

To improve upon the prediction of Eq. ([]]) for the region Q T <C Q, perturbation theory 
can be applied using an expansion parameter a™ In" (Q 2 j Q T ) , with n — 0, . . . ,2m — 1, 
instead of a™. The terms a™ ln n (Q 2 /Q' T ) represent the effects of soft gluon emission 
at order a™. Resummation of the singular part of the perturbative series to all orders 
in a s by Sudakov exponentiation yields a regular differential cross section as Qt — > 0. 

The differential cross section in the CSS resummation formalism for the production 
of photon pairs in hadron collisions can be written as an integral in impact parameter 
b space: 

do(hih 2 — ► 7172^0 1 1 

dQ 2 dydQ 2 T d cos 9 d<p ~ 24nS Q 2 

1- J d 2 be i ^ rb Y,W t3 (K,Q,x 1 ,x 2 ,6,<p,C u C 2 ,C 3 )W^ p (b,Q,x 1 ,x 2 ) 



x 



1 -j 



+ Y(Q T ,Q,x u x 2 ,6,<j ) ,C A )}. (2) 

The variables Q, y, and Qt are the invariant mass, rapidity, and transverse momentum 
of the photon pair in the laboratory frame, while 9 and <ft are the polar and azimuthal 



angle of one of the photons in the Collins-Soper frame [ 18| . The initial-state parton 
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momentum fractions are defined as x\ = e y Q/\/~S, and x 2 = e~ y Q/\fS, and \fS is 
the center-of-mass (CM) energy of the hadrons h\ and h 2 . 

The renormalization group invariant quantity Wij(b) sums the large logarithmic 
terms a™ln n (6 2 Q 2 ) to all orders in a s . For a hard scattering process initiated by the 
partons i and j, 

Wij(b, Q, xi, x 2 , 9, 4>, Ci, C 2 , C 3 ) = exp {-<Sy(&, Q, Ci, C 2 )} 



x 



Ci/ hl (xx)Cj/h 2 {x 2 ) +C j / hl (x 1 )Ci/ h2 (x 2 ) Fij(a em (C 2 Q) , a s (C 2 Q) , 9, 0), (3) 



where jFy is a kinematic factor that depends also on the coupling constants, and 
Ci/h(x) denotes the convolution of the parton distribution function (PDF) f a / h (for 
parton a inside hadron h) with the perturbative Wilson coefficient functions C^: 

C i/hl (*i) = £ £ f 1 (|, M = y , G.) /./* A* = ■ (4) 

The Sudakov exponent Sij(b,Q,Ci,C 2 ) in Eq. (§) is defined as 

c|Q 2 ci/i 2 



S tJ (b,Q,C u C 2 )= f 



Cf/b 2 /i 2 



Aij (a,(fi), Ci) In + % (a s (/x), C 1; C 2 ) 



(5) 

The coefficients Ay and .By and the functions Cy can be calculated perturbatively 

oo oo 

in powers of a s /ir, so that Ay = £(a s /7r) n A^; , -By = y^(a s /7r) n i?y"\ and Cy = 



n=l n=l 

t(n) 



E(«./t) b C- 

n=0 

The dimensionless constants Ci, C 2 and C3 = /i& were introduced in the solution 
of the renormalization group equations for Wij. The constant C\ determines the onset 
of nonperturbative physics, C 2 specifies the scale of the hard scattering process, and 
jj, = C?,/b is the factorization scale at which the functions are evaluated. A 
conventional choice of the renormalization constants is C\ — C 3 — 2e~ lE = b and 
C 2 = C4 = 1 [[/]], where 7^ is the Euler constant. These choices of the renormalization 
constants are used in the numerical results of this work because they eliminate large 
constant factors (depending on Ci,C 2 and C 3 ) in the Sudakov exponent and in the 
C\f functions 0. 
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In Eq. (fj), the impact parameter b is to be integrated from to oo. However, for 
b > b max , which corresponds to an energy scale less than l/5 max , the QCD coupling 
a s (jl ~ 1/b) becomes so large that a perturbative calculation is no longer reliable, and 
nonperturbative physics must set in. The nonperturbative physics in this region is 
described by an empirically fit function W^ p || ^0| , and is evaluated at a revised 
value of b, 6* = , where 6 max is a phenomenological parameter used 

to separate long and short distance physics. With this change of variable, 6* never 
exceeds 6 max ; b max is a free parameter of the formalism M that can be constrained by 
other data (e.g. Drell-Yan). 

The function Y in Eq. (0) contains contributions in the full NLO perturbative 
calculation that are less singular than Q^ 2 or Q^ 2 \n(Q 2 /Q T ) as Qt — > (both 
the factorization and the renormalization scales are chosen to be C±Q). It is the 
difference between the exact perturbative result to a given order and the result from 
Wij expanded to the same fixed order (called the asymptotic piece). The function 
Y restores the regular contribution in the fixed order perturbative calculation that 
is not included in the resummed piece W^. It does not contain a contribution from 
final-state fragmentation, which is included separately as described in Sec. 3.2. 

The CSS formula Eq. (0) contains many higher-order logarithmic terms, such 
that when Qt ~ Q, the resummed differential cross section can become negative in 
some regions of phase space. In this calculation, the fixed-order prediction for the 
differential cross section is used for Qt ~ Q whenever it is larger than the prediction 
from Eq. (|j). The detailed properties of this matching prescription can be found in 
Ref. 0. 

3.1 Resummation for the qq —> 77 subprocess 

For the qq —>■ 77 subprocess, the application of the CSS resummation formalism 
is similar to the Drell-Yan case qq ^ — > V* — > £^2 , where i\ and £2 are leptons 
produced through a gauge boson V* |10| . Since both processes are initiated by qq ^ 



color singlet states, the A^\A^ 2 > and functions in the Sudakov form factor are 



7 



identical to those of the Drell-Yan case when each photon is in the central rapidity 
region with large transverse momentum and is well separated from the other photon. 
This universality can be understood as follows. The invariants s, i and u are defined 
for the q{pi)q{p-i) — > lip^lipi) subprocess as 

s = (pi +p 2 ) 2 , t = (p 1 -p 3 ) 2 , u=(p 2 -p 3 ) 2 . (6) 



The transverse momentum of each photon can be written as p T = y tuj s. When p\ is 
large, t and u must also be large, so the virtual-quark line connecting the two photons 
is far off the mass shell, and the leading logarithms due to soft gluon emission beyond 
the leading order can be generated only from the diagrams in which soft gluons are 
connected to the incoming (anti-)quark. To obtain the B^ function, it is necessary 
to calculate beyond NLO, so it is not included in this calculation. However, the 
Sudakov form factor becomes more accurate when more terms are included in 
and Bij. Since the universal functions A^" depend only on the flavor of the incoming 
partons (quarks or gluons), A q q can be appropriated from Drell-Yan studies, and its 
contribution is included in this paper. 

To describe the effects of multiple soft-gluon emission, Eq. (|2|) can be applied, 
where i and j represent quark and anti-quark flavors, respectively, and Tij = 5ij(g\ + 
<7#) 2 (1 + cos 2 0)/{l — cos 2 6). The couplings gi,n are defined through the qqj vertex, 
written as [g L (l - 75) + ^(1 + 75)] , with g L = g R = eQ f /2, and eQ f is the 
electric charge of the incoming quark (Q u = 2/3, Qd = —1/3). The explicit forms of 
the A and B functions are: 



Aqq (Ci) — Cp, 



F 







Bq\ (Cl; C2) — Ci 



.1-3 

3 01 (C 2 b» 
— I In 



2 -y Cl • ■ (7) 

where Nf is the number of light quark flavors, Nc = 3, Cp = 4/3, and (3\ = (UNc — 
2JV»/12. 

To obtain the value of the total cross section to NLO, it is necessary to include 
the Wilson coefficients Cl° and C\)\ These can be derived from the full set of LO 



contributions and NLO corrections to 77 production [H. After the leading order and 
the one-loop virtual corrections to qq — > 77 and the tree level contribution from 
99 — * 119 are included, the coefficients are: 

C$(z,b,fi,-jj-) = 6 jk 5(l-z), 



0. 



+5(1-2) 



& C 2 



V _9 

4 16 



After factorization of the final-state collinear singularity, as described below, the real 
emission subprocess qg — »■ 775 yields: 



tv2 



In the above expressions, the splitting kernels [19 



are 



1 + z 2 
1 - z 



z 2 + (l-z) 



and 



(9) 



(10) 



For photon pair production, the function V is 

r 2 



V 



77 



~ 4 + T + fl^TP ( F °"*( v ) - 2 



ln 2 (w) 



+ 



1 - v 
V 1 — V 



1 - V 



ln 2 (l - v) 



l-v 



V 



\n 2 (v) + ln 2 (l - v) - 3 



2vr 2 



+2 (ln(u) + ln(l - u) + 1) + 3 ( — — ln(l - v) + - — - ln(v 

VI — v v 



(11) 



where v = —u/s, and u = —s(l + cos#)/2 in the gg center-of-mass frame. Because 
of Bose symmetry, F mrt {v) = F mrt (l — v). A major difference from the Drell-Yan 
case (Vdy = — 8 + rr 2 ) is that V 77 depends on the kinematic correlation between the 
initial and final states through its dependence on u and t. 
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The nonperturbative function used in this study is the empirical fit p][ 



W"«Sj P ( 6 > Q> Qo> x ^ x i) = ex P 



-gib 2 - g 2 b 2 In ( J - gxg 3 bla (100xix 2 ) 



(12) 



2Q , 

where g x = 0.11±g;g| GeV 2 , g 2 = 0.581^ GeV 2 , g 3 = — 1.5±g;l GeV~\ and Q = 
1.6 GeV. (The value 6 max = 0.5 GeV -1 was used in determining the above g^s 
and for the numerical results presented in this paper.) These values were fit for the 
CTEQ2M parton distribution function, with the conventional choice of the renormal- 
ization constants, i.e. C\ = C3 = bo and C 2 = 1. In principle, these coefficients 



should be refit for the CTEQ4M distributions |2I] used in this study. The parameters 
of Eq. flT2"D were determined from Drell-Yan data. It is assumed that the same values 
should be applicable for the 77 final state. 

3.2 Contributions From qg Subprocesses 

As described in Sec. 2, the complete NLO calculation of diphoton production in 
hadron collisions includes photons from long-distance fragmentation processes like 
(lb) and short-distance processes like (If) and (lg). The latter processes yield a 
regular 3-body final state contribution, while the former describes a photon recoiling 
against a collinear quark and photon. 

The singular part of the squared amplitude of the q(pi)g(p 2 ) — > 7(^3)7(^4)9(^5) 
subprocess can be factored into a product of the squared amplitude of q(j>i)g(p 2 ) —* 
7(^3)9(^4+5) an d the splitting kernel for 5(^4+5) — > 7(^4)9(^5) • In the limit that the 
emitted photon 7(^4) is collinear with the final state quark q{jp§): 

Km \M (q(pi)g(p 2 ) -> 7(^3)7(^4)9(^5)) I 2 = 



,2 

— I 

7<-q 



P« (z) \M (q(pi)g(p 2 ) - 7(P3)9(P4 +5 ))| 2 • (13) 



Pi -Pb 

A similar result holds when p% and p§ become collinear and/or the quark is re- 
placed with an anti-quark. Conventionally, the splitting variable z is the light-cone 
momentum fraction of the emitted photon with respect to the fragmenting quark, 
z — P4 / {pi + vt)i where pf = [pf^ + pf'^j /y/2. (p^ is the energy and pf^ is the 
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longitudinal momentum component along the moving direction of the fragmenting 
quark in the qg center-of-mass frame.) Alternatively, since the final state under con- 
sideration contains only a fragmenting quark and a spectator, a Lorentz invariant 



splitting variable can be defined as: [21] 



~z = \ . (14) 

Pj -Pk+Pi- Pk + Pi ■ Pj 

In this notation, i = 5 is the fragmentation quark, j is the fragmentation photon, and 
k is the prompt spectator photon. When the pair ij becomes collinear, z becomes 
the same as the light-cone momentum fraction z carried by the photon. Aside from 
the color factor, P^}_ q {z) in Eq. (|1|) is the usual DGLAP splitting kernel for q — > gq 

po, w = (UllzfJI) . (15) 

The regular contribution qg —>■ 775 (If) is defined by removing the final-state, 
collinear singularity from the full amplitude of the partonic subprocess. The matrix 
element squared for (If) can be written |[21|| : 



\M {qg -» 77?) | 



2 

reg 



\M {qg - liq )\) ull - (*) \M (g(Pi)y(Pa) - i(pMP4 + ,))\ 2 ■ (16) 

After the final-state collinear singularity is subtracted, the remainder expresses the 
regular 3-body final state contribution 77?. This remainder, as shown in (lg), con- 
tains terms that diverge when Qt — > which should be regulated by renormalizing 
the parton distribution at the NLO. The contribution from this divergent part is in- 
cluded in the resummed qq cross section in Cjq, as shown in Eq. (|]). The part that 
is finite as Qt — > is included in the function Y. When Qt ~ Q, Eq. (p~6|) describes 
the NLO contribution from the qg — ► 77? subprocess to the Qt distribution of the 
photon pair. The subtracted final-state collinear singularity from the NLO qg — ► 77? 
subprocess is absorbed into the fragmentation process (lb). 

3.3 Fragmentation Contributions 

Final-state photon fragmentation functions D 7 /i(z, fi 2 F ) are introduced in an analo- 
gous manner to initial-state parton distribution functions fi/hifai lA)- Here, z{x) is 
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the light-cone momentum fraction of the fragmenting quark (incident hadron) carried 
by the photon (initial-state parton), and is the final state (initial state) frag- 

mentation (factorization) scale. The parton-level cross section for the fragmentation 
contribution (lb) is evaluated from the general expression for a hard scattering to a 
parton m, which then fragments to a photon: 

da=^\M(p lP2 ^p 3 ...p m )\ 2 d^[PS]dzD y/m (z,fi 2 F ). (17) 

Here, Ai is the matrix element for the hard scattering subprocess, S m ^[PS] is the 
m-body phase space, and an integral is performed over the photon momentum frac- 
tion z weighted by the fragmentation function D 7 / m (z, ^ F ). Since fragmentation is 
computed here to LO only, the infrared divergences discussed by Berger, Guo and 



Qiu are not an issue [22 



The fragmentation function D^^ q obeys an evolution equation, and the leading- 
logarithm, asymptotic solution D^f q is 0: 



QCD 

Q 2 J2.21-1.28z + 1.29z 2 )z 0M9 , x9n 1 ,, . 

zD ^ = 8 i- 1.63m-,) +°- roa '( 1 -»)"»-"•■ < lf 



where Aqcd is the QCD scale for four light quark flavors. As shown in Fig. 0, the 
collinear approximation made in defining D^^ q leads to kinematic distributions with 
an unrealistic sensitivity to kinematic cuts, such as cuts to define an isolated photon. 

The Monte Carlo showering method goes beyond the collinear approximation used 
in solving the evolution equation for the fragmentation function -D 7 <_ g . In Monte Carlo 
calculations, the probability for photon emission is determined from the splitting 
function P y ^ q (z), which is a collinear approximation. However, the kinematics are 
treated by assigning a virtuality to the fragmenting quark whose value lies between 
the hard scale of the process and a phenomenological cutoff ~ 1 GeV. This cutoff 
replaces the parameter Aqcd in Eq. flTE|). Most importantly, gluon emission can 
be incorporated into the description of final state fragmentation. Because there is 
no collinear approximation in the kinematics, kinematic distributions do not exhibit 
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the unrealistic behavior of the parton-level calculation. The "correctness" of either 
approach can be judged only after a careful comparison of their respective predictions. 

The collinear approximation becomes an issue because of the experimental def- 
inition of isolated photons. Experimentally, an isolation cut is necessary to sepa- 
rate prompt photons from various hadronic backgrounds, including tt° and rj me- 
son decays. The separation between a particle j and the photon is expressed as 



= V ~~ ^i) 2 + ~ ^i) 2 ' where the coordinates r)(rjj) and 4>{4>j) are the pseu- 
dorapidity and azimuthal angle of the photon (particle j). At hadron colliders, the 
standard isolation criterion is that the sum of excess transverse energy Et contained 
inside a cone of size Rq centered on the photon candidate is below a cutoff E^°, 
E T < E^°. The sum is over each particle j. Since the resummed CSS piece 

Rj<Ro 

of the final state cross section describes the radiation of multiple soft gluons approx- 
imately collinear with the incident partons, it produces only isolated photons. For 
NLO 77 j final states (le), (If), and (lg), where there is only one extra parton j = q 
or g, isolation enforces a separation Rj > Rq, provided that Pt > E^°. Above 
Qt = E™ , the perturbative corrections contained in the function Y are affected by 
isolation. On the other hand, because of the collinear approximation, the parton- 



level fragmentation calculation based on Eq. (18) does not depend on the isolation 
cone Rq; the hadronic remnant of the fragmentation (lb) always satisfies R < Rq. 
Hence, for this case, Qt = (1 — z)pr 1 , and the isolation cut reduces to a step function 
requirement 9(E^° — Qt)- 

The parton-level calculation of the fragmentation contribution at the Tevatron 
based on the fragmentation function D^^ q (z, //f,) has been compared with a Monte 



Carlo estimate based on PYTHIA [p3| . For the parton-level calculation, the scale ftp = 
M 77 is used. For the PYTHIA calculation, the scale is pip — y/I, and hadronization is 
not performed, so that no photons arise from 7r° or rj meson decays, for example. For 
this comparison, the invariant mass of the hard-scattering subprocess is limited to 
20 < \/§ < 50 GeV in both approaches, and the photons are required to satisfy p T > 5 
GeV and |?] 7 | < 2. These kinematic cuts are chosen to increase the statistics of the 
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PYTHIA calculation, while reflecting the kinematic region of interest for a comparison 
with data. PYTHIA can simulate the QED and QCD showering of the final-state quark 
as well as the QCD showering of the initial-state quark and gluon. To isolate the effect 
of initial-state gluon radiation, PYTHIA calculations were performed with and without 
the QCD initial-state radiation (i.e. by preventing space-like showering). In neither 
case is initial-state QED radiation simulated. It is possible for the partons produced 
in initial-state showering to develop time-like showering. Any photons produced 
from this mechanism are discarded, since they are formally of higher-order than the 
contributions considered here. Such contributions, however, might be necessary to 
understand photon pairs with small invariant mass and small Qt- 

Figure Q is a comparison of kinematic quantities from the parton-level and Monte 
Carlo calculations. The left-side of Fig. ||] shows the (^-distribution for the parton- 
level (solid), PYTHIA with initial-state radiation of gluons (short-dashed), and PYTHIA 
without initial-state radiation (long-dashed) calculations. Each curve is plotted 
twice, with and without an isolation cut = 4 GeV and Rq = 0.7. Before the 
isolation cut, the total parton-level fragmentation cross section is approximately 50% 
higher than the Monte Carlo cross section. After isolation, the total cross sections 
are in good agreement, even though the parton-level calculation is discontinuous at 
Qt = E™ . The effect of initial-state gluon radiation in the PYTHIA calculation is to 
increase Qt without compromising the isolation of the photons. 

The right-side of Fig. Q shows the distribution of the light-cone momentum frac- 
tion z of the quark carried by the fragmentation photon (for this figure, z is defined 
in the laboratory frame). After isolation, the parton-level contribution is limited to 
z > 0.55 by kinematics, whereas the Monte Carlo contribution is more uniformly dis- 
tributed between and 1. For the PYTHIA result, z is calculated with respect to the 
final state quark before showering. In the showering process, some energy-momentum 
can be exchanged between the final state prompt photon and the fragmenting quark, 
since the quark is assigned a virtuality. As a result, the effective z value can extend 
beyond the naive limit z = 1. 
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The conclusions of this comparison are: (1) after isolation, the total cross sec- 
tions from the parton-level and Monte Carlo fragmentation calculations are in good 
agreement, and (2) the Monte Carlo kinematic distributions (e.g. Qt and z) are not 
very sensitive to the isolation cut. For these reasons, the Monte Carlo estimate with 
initial-state radiation is used to account for the (lb) contribution in the final results. 
Furthermore, with initial-state radiation, the PYTHIA calculation includes the leading 
effects of a full resummation calculation of the qg — > jq process. It is approximately 
equivalent to performing a resummation calculation in the CSS formalism with quan- 
tities and B^> calculated for a qg initial state and the LO Wilson function. 

One final comparison was made with the Monte Carlo calculation by treating 
the subtracted term in Eq. (16), with replaced by defined in Eq. (18), 
as a 3-body matrix element. The collinear divergence was regulated by requiring a 
separation Rq between the photon and quark remnant for all Qt- This calculation 
agrees with PYTHIA in the shape and normalization of various distributions, except 
when Q T < E™ , where there is a substantial difference. 

3.4 Resummation for the gg — > 77 subprocess 

A resummation calculation for the gg — > 77 subprocess is included in the theoretical 
prediction. The LO contribution comes from one-loop box diagrams of order a 2 em a 2 s in 
perturbative QCD. At present, a full NLO calculation, of 0(al m a^), for this process is 
not available. Nevertheless, the resummation technique can be applied to resum part 
of the higher order contributions and improve the theoretical prediction. The exact 
NLO gg — > 77$ calculation must include gluon emission from the internal quark lines 
of the box diagram, thus generating pentagon diagrams. However, such diagrams 
do not generate large logarithms when the final state photons have large transverse 
momentum, are in the central rapidity region, and are well separated from each other. 
All the large logarithms originate from the diagrams with soft gluons coupling to the 
initial-state gluons. Similarly, the exact NLO qg — > 779 calculation, of 0(al m a^), 
must include contributions involving a box diagram with one incoming gluon off 
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shell. Large logarithms only arise from soft gluon emission off the initial-state quark 
or gluon. The leading logarithms due to initial-state radiation are universal, and the 
function calculated for the resummed gg — > H process |fLl , or the color singlet 



part of the gg — > QQ process can be applied directly to the resummed gg — > 77 
calculation, since these subprocesses have the same QCD color structure. 

When the transverse momentum of the photon pair is much smaller than its 
invariant mass, i.e. Qt <C Q, and each photon has large transverse momentum, then 
the box diagram of the hard scattering subprocess gg — *> 77 can be approximated 
as a point-like interaction (multiplied by a form factor which depends on s, t and 
u). This approximation ignores pentagon diagrams in the gg — > 77*7 subprocess and 
the virtuality of intermediate quarks in the qg — > 77^ subprocess. It does not have 
the complete structure of the hard process, but it does contain the most important 
logarithmic terms from initial state gluon radiation. Under such an approximation, 
the subleading logarithmic terms associated with A^ 2 \ and of Eqs. (4) and 
(5) can be included in the resummation calculation. These functions were calculated 
for the gg — > H process JTTJ Without a complete 0(al m a^) calculation, the 
exact Wilson coefficient function is not known. Since part of the exact 
function must include the piece for the gg — > H process, it is included to estimate the 
possible NLO enhancement to the production rate of the gg subprocess. After these 
ingredients are incorporated into Eq. (Q), the resummed kinematics of the photon 
pair from the gg — > 77 subprocess can be obtained. The distribution of the individual 
photons can be calculated approximately from the LO angular dependence of the box 
diagram. 

The above approximation certainly fails when Qt is of the order of Q. In the 
absence of a complete 0{a1 m a z s ) calculation of the gg — > 77$ and qg — > 779 subpro- 
cesses, it is not possible to estimate the uncertainties introduced by the approxima- 
tion. In the limit of Qt <^ Q, the approximation should be reliable, since the soft 
gluon approximation is applicable. In the same spirit, the approximate function Y 
for photon pair production is taken from the results of the perturbative piece for the 



16 



gg — > Hg and gq — > Hq processes [|TT], |T^]. 

In summary, the resummed distributions of the photon pair from the gg subpro- 
cess in the region of Qt <C Q can be described by Eq. (|2|), with i — j — g, and 
T gg = N c \M.gg^ n {s,t,u)\^ /2 12 . Here, \M. gg ^ T( {s, t, u)\ 2 is the absolute square of 
the invariant amplitude of the gg — > 77 subprocess || summed over spins, colors, 
and the fermion flavors in the box loop, but without the initial-state color (1/8 2 ), 
spin (1/2 2 ) average, and the final-state identical particle (1/2) factors. The A and B 
functions used in the calculation for the gg initial state are 
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H subprocess, are: 
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(20) 



Since the NLO pentagon and off-shell box diagram calculations are not included, 
the Wilson coefficients C\f are expected to predict accurately the total cross section 
only when Qt -C Q, the transverse momenta of the individual photons are large, and 
their rapidities are small. Under the approximation made above, the resummed gg 
result increases the integrated rate by about a factor of 2, for kinematic cuts typical 
of the Tevatron, as compared to the lowest order (one-loop calculation) perturbative 
result. This comparison suggests that the full NLO contribution of the gg initiated 
subprocess is large. Because it is necessary to impose the condition Qt < Q to 
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make the above approximations valid, the gg resummed result presented in this work 
probably underestimates the rate when Qt is large or the separation of the azimuthal 
angle (A0) between the two photons is small. This deficiency can be improved only 
by a complete 0(a^ m a^) calculation. 

At the Tevatron, the gg contribution is important when the invariant mass (M 77 = 
Q) of the two photon pair is small. Because of the approximation made in the gg 
calculation beyond the LO, the prediction will be more reliable for the data with 
larger Q. A more detailed discussion is presented in the next section. 

The full calculation of the gg contribution in the CSS formalism depends also upon 
the choice of nonperturbative functions. However, the best fits to the parametriza- 
tions are performed for qq processes P, EDI. Two assumptions were studied: (i) the 



nonperturbative functions are truly universal for qq and gg processes, and (ii) the 
nonperturbative functions obey the same renormalization group properties as the A 
functions for each type of process (which are universal for all qq or gg subprocesses), so 
the coefficient of the In (2^) term in the nonperturbative function Eq. (|l^) is scaled 
by Ca/Cf relative to that of the qq process. Specifically, the different assumptions 
are: 

(i)W g N g p (b,Q,Qo,x 1 ,x 2 ) = W q p i p (b,Q,Q ,x 1 ,x 2 ) 



(ii)W g N g p (b,Q,Q ,x 1 ,x 2 ) = W^ g p (b,Q,Q ,x 1 ,x 2 )(g 2 ^^g 2 ) 



(21) 



The numerical values of gi,g 2 , and #3 are listed following Eq. (12). These two as- 
sumptions do not exhaust all possibilities, but ought to be representative of reason- 
able choices. Choice (ii) is used for the results presented in this paper. The effect of 
different choices is discussed in Sec. 4. 

4 Numerical Results 
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4.1 Tevatron Collider Energies 

Two experimental collaborations at the Tevatron pp collider have collected diphoton 
data at = 1.8 TeV: CDF fl|], with 84 ptr\ and D0 fTJ, with 81 pb _1 . The 
kinematic cuts applied to the resummed prediction for comparison with the CDF data 
are p\ > 12 GeV and |^ 7 | < 0.9. For D0, the kinematic cuts are p^ > 14 GeV and 
Pt? > 13 GeV, and |7y 7 | < 1. For CDF, an isolation cut for each photon of R = 0.7 
and = 4 GeV is applied; for D0, the cut is R = 0.4 and E l T so = 2 GeV. 

Other ingredients of the calculation are: (i) the CTEQ4M parton distribution func- 
tions, (ii) the NLO expression for a s , (in) the NLO expression for a em , and (iv) the 
nonperturbative coefficients of Ladinsky-Yuan 



The predictions for the CDF cuts and a comparison to the data are shown in 
Figs. |3|-|5]. Figure [3] shows the distribution of the photon pair invariant mass, da/dM^ 
vs. M 77 . The dashed-dot curve represents the resummation of the gg subprocess, 
which is the largest contribution for M 77 ~ 30 GeV. The long-dashed curve represents 
the full qq resummation, while the short-dashed curve is a similar calculation with 
the gluon parton distribution function artificially set to zero. Schematically, there are 
contributions to the resummed calculation that behave like q — > gqi <S> qiq — > 77 and 
g — > qqi g) q\q — > 77. These contributions are contained in the terms proportional 
to Pfl k (z) in Eq. (§) and Pfl G (z) in Eq. (|D|), respectively. The full qq resummation 
contains both the qq and qg contributions. The short-dashed curve is calculated by 
setting Cjq = and retaining only the qq contribution in the function Y. Since the 
short-dashed curve almost saturates the full qq + qg contribution, except at large 
Qt or small A<f>, the qg initiated subprocess is not important at the Tevatron in 
most of phase space for the cuts used. The fragmentation contribution is denoted 
by the dotted line. The sum of all contributions including fragmentation is denoted 
by the solid line. After isolation, the fragmentation contribution is much smaller 
than "direct" ones, but contributes ~ 10% near the peak. The uncertainty in the 
contribution of the fragmentation process can be estimated by comparing the Monte 
Carlo result with a parton-level calculation, as shown in Fig. |[ 
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Figure |] shows the distribution of the transverse momentum of the photon pair, 
da/dQx vs. Qt- Over the interval 5 ~ Qt ~ 25 GeV, the contribution from the gg 
subprocess is comparable to the qq+qg subprocess. The change in slope near Qt = 20 
GeV arises from the gg subprocess (dot-dashed line) for which Qt ~ M 77 is required 
in our approximate calculation. The peak near Qt ~ 1.5 GeV is provided mostly by 
the qq + qg (long-dashed line) subprocess. In general, the height and breadth of the 
peak in the Qt distribution depends on the details of the nonperturbative function in 
Eq. (0). The effect of different nonperturbative contributions may be estimated if the 
parameter g 2 is varied by ±2cr. For Qt > 10 GeV, the distribution is not sensitive 
to this variation. The height and the width (half-maximum) of the peak change 
by approximately 20% and 35%, respectively, but the integrated rate from to 10 
GeV is almost constant. The peak of the distribution (which is below 5 GeV), shifts 
approximately +0.5 GeV or —0.6 GeV for a +2<r or —2a variation. The mean Qt for 
Qt < 10 GeV shifts at most by 0.4 GeV. For gg resummation, it is not clear which 
parametrization of the nonperturbative physics should be used. However, the final 
effect of the two different parametrizations outlined in Eq. fl2T|) is minimal, shifting 
the mean Qt for Qt < 40 GeV by about 2.0 GeV. The parametrization (ii) is used 
in the final results, so that the coefficient g 2 is scaled by Ca/Cf relative to the qq 
nonperturbative function. 

Figure |5| shows da/dA<p vs. A<p, where A0 is the azimuthal opening angle between 
the two photons. The change in slope near A0 = it/ 2 is another manifestation of 
the approximations made in the treatment of the gg contribution (dot-dashed line). 
The height of the distribution near A0 ~ tt is also sensitive to the details of the 
nonperturbative function. 

In the absence of resummation or NLO effects, the gg box contribution supplies 
Qt = and A0 = tt. In this calculation, as explained earlier, the NLO contribution 
for the gg subprocess is handled in an approximate fashion. For the cuts listed above, 
the total cross section from the complete gg resummed calculation, including the 
function Y, is 6.28 pb. If the resummed CSS piece is used alone, the contribution is 
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reduced to 4.73 pb. This answer can be compared to the contribution at LO. For the 
same structure functions, the LO gg cross section for the CDF cuts is 3.18 pb for the 
scale choice Q = M 77 . Therefore, the effect of including part of the NLO contribution 
to the 77 process is to approximately double the LO gg box contribution to the cross 
section. This increase indicates that the exact NLO correction can be large for the 
gg subprocess and motivates a full calculation. 

The predictions for the D0 cuts and a comparison to data are shown in Figs. 
[|. Because of the steep distribution in the transverse momentum of the individual 
photons, the higher p T threshold in the D0 case significantly reduces the total cross 
section. Otherwise, the behavior is similar to the resummed calculation shown for 
the CDF cuts. The D0 data plotted in the figures are not corrected for experimental 
resolution. To compare with the uncorrected D0 data with the kinematic cuts p T x > 
14GeV, pj 2 > 13GeV and rf < 1.0, an "equivalent" set of cuts is used in the 



theoretical calculation: p T x > 14.9 GeV, p£ > 13.85 GeV, and if < 1.0 [||. The 
effect of this "equivalent" set is to reduce the theoretical rate in the small M 77 region. 

While the agreement in both shapes and absolute rates is generally good, there 
are some discrepancies between the resummed prediction and the data as presented in 
these plots. At small Qt (Fig. |j) and large (Fig. where the CDF cross section 
is large, the theoretical results are beneath the data. Since this is the kinematic region 
in which the nonperturbative physics is important, better agreement can be obtained 
if the nonperturbative function is altered. In Fig. |^, the calculated M 77 distribution 
is larger than the D0 data at large M 77 , while the calculation appears to agree with 
the CDF data in Fig. [|[ The small discrepancy in Fig. ^ at large values of M 77 is 
not understood. (The systematic errors of the data, which are about 25% f25| , are 



not included in this plot.) On the other hand, Figs. [7| and [8] show that the resummed 
calculation is beneath the data at large Qt or small A<fr. The discrepancies in Figs. [F] 
and H] may result from the approximations made in the gg process (notice the kinks in 
the dot-dashed curves). A complete NLO calculation for the gg subprocess is needed, 
and may improve the comparison with data for small A0. 
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Because of the uncertainty in the prediction for the gg contribution of the re- 
summed calculation, the distributions in M 77 , Qt and A0 are shown in Figs. |D] |ll 



for the CDF cuts and the additional requirement that Qt < M 77 . This additional 
requirement should significantly reduce the theoretical uncertainty for large Qt and 
small A(j). 

In Fig. [K| the lower of the two solid curves in the Qt distribution shows the 
prediction of the pure NLO 0(a s ) (fixed-order) calculation, without resummation, 
for the qq and qg subprocesses, excluding fragmentation. For Qt ~ 25 GeV, the lower 
solid curve is very close to the long-dashed (qq + qg) curve obtained after resumma- 
tion, as is expected. As Qt decreases below Qt — 25 GeV, all-orders resummation 
produces significant changes. Most apparent, perhaps, is that the Qt — > divergence 
in the fixed-order calculation is removed. However, there is also a marked difference 
in shape over the interval 5 < Qt < 25 GeV between the fixed-order qq + qg re- 
sult and its resummed counterpart. These are general features in a comparison of 



resummed and NLO calculations HI— [14 



4.2 Fixed— Target Energy 



The fixed-target experiment E706 |17j at Fermilab has collected diphoton data from 
the collision of a p beam on a Be (A = 9.01, Z = 4) target at v 7 ^ = 31.5 GeV. The 
kinematic cuts applied to the resummed prediction in the center-of-mass frame of 
the beam and target are p T > 3 GeV and |^ 7 | < 0.75. No photon isolation is required. 
The same phenomenological inputs are used for this calculation as for the calculation 
at collider energies. The Be nucleon target is treated as having an admixture of 4/9.01 
proton and 5.01/9.01 neutron parton distribution functions. The A dependence effect 
appears to be small in the prompt photon data (the effect is parametrized as A a and 
the measured dependence is a ~ 1), so it is ignored P5| . 

Figures |l2HT4l show the same distributions discussed previously. Because of the 
kinematic cuts, the relative contribution of gluon initiated processes is highly sup- 
pressed, except at large Qt, where the gg box contribution is seen to dominate, and 
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at large M 77 where the qg contribution is dominant. The fragmentation contribution 
(not shown) is minimal (of a few percent). The dominance of gg resummation over 
the qq resummation at large Qt in Fig. [13] occurs because it is more likely (enhanced 
by the ratio Ca/Cf = 9/4) for a gluon to be radiated from a gluon than a quark line. 
The exact height of the distribution is sensitive to the form of the nonperturbative 
function (in the low Qt region) and to the approximation made in calculating the 
NLO corrections (of 0(al m a^)) to the hard scattering. However, since Qt < Q is 
satisfied for the set of kinematic cuts, the final answer with complete NLO corrections 
should not differ significantly from the result reported here. 

The scale dependence of the calculation was checked by comparing with the result 
obtained with C 2 = Ci/b = 0.5, C 3 = b , and C 4 = 1. The qq rate is not sensitive 
to the scale choice, and the gg rate increases by less than about 20%. This insensi- 
tivity can be understood as follows. For the E706 data, the nonperturbative physics 
completely dominates the Qt distribution. The perturbative Sudakov resummation 
is not important over the entire Qt region, and the NLO Y piece is sizable only for 
Qt > 3 GeV where the event rate is small. Since the LO qq rate does depend on a s , 
and the LO gg rate is proportional to a^(C2M 77 ), the gg rate increases for a smaller 
C*2 value, but the qq rate remains about the same. In conclusion, the E706 data can 
be used to constrain the nonperturbative functions associated with the qq and gg 
hard processes in hadron collisions. 

5 Discussion and Conclusions 

Prompt photon pair production at fixed target and collider energies is of interest in its 
own right as a means of probing the dynamics of strong interactions. The process is 
of substantial interest also in searches for new phenomena, notably the Higgs boson. 

In this paper, a calculation is presented of the production rate and kinematic 
distributions of photon pairs in hadronic collisions. This calculation incorporates the 
full content of the next-to-leading order (NLO) contributions from the qq and qg 
initial-state subprocesses, supplemented by resummation of contributions to these 
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subprocesses from initial state radiation of soft gluons to all orders in the strong 
coupling strength. The computation also includes important contributions from the 
gg box diagram. The gg contributions from initial-state multiple soft gluons are re- 
summed to all orders, but the NLO contribution, of 0{a 2 em a z s ), to the hard scattering 
subprocess is handled in an approximate fashion. The approximation should be reli- 
able at relatively small values of the pair transverse momentum Qt as compared to 
the invariant mass of the photon pair M 77 . At collider energies, the gg contribution 
is comparable to that of the qq and qg contributions over a significant part of phase 
space where M 77 is not large, and its inclusion is essential. The exact 0(al m al) 
corrections to the gg box diagram should be calculated to test the validity of the 
approximations made in this calculation. Finally, the calculation also includes long- 
distance fragmentation contributions at leading order from the subprocess qg — > jq, 
followed by fragmentation of the final quark, q — > jX. After photon isolation, frag- 
mentation plays a relatively minor role. The fragmentation contribution is computed 
in two ways: first, in the standard parton model collinear approximation and second, 
with a Monte Carlo shower simulation. This overall calculation is the most complete 
treatment to date of photon pair production in hadronic collisions. Resummation 
plays a very important role particularly in the description of the behavior of the Qt 
distribution at small to moderate values of this variable, where the cross section takes 
on its largest values. 

The resummed calculation is necessary for a reliable prediction of kinematic distri- 
butions that depend on correlations between the photons. It is a significant improve- 
ment over fixed-order NLO calculations that do not include the effects of initial-state 
multiple soft-gluon radiation. Furthermore, even though the hard scattering qq and 
qg subprocesses are computed to the same order in the resummed and fixed-order 
NLO calculations, the cross sections from the two calculations can differ after kine- 
matic cuts are imposed [ |10|| . 

The results of the calculation are compared with data from the CDF and D0 col- 
laborations, and the agreement is generally good in both absolute normalization and 
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shapes of the distributions in the invariant mass M 77 of the diphoton system, the pair 
transverse momentum Qt, and the difference in the azimuthal angles A0. Discrep- 
ancies with CDF results at the smallest values of Qt and A</> near it might originate 
from the strong dependence on the nonperturbative functions in this kinematic re- 
gion. In comparison with the D0 data, there is also evidence for disagreement at 
intermediate and small values of A0. The region of intermediate A0, where the two 
photons are not in a back-to-back configuration, is one in which the full treatment 
of three body final-state contributions of the type 77 j are important, with j = q or 
g. The distributions in Figs. 5 and 8 suggest that an exact calculation of the NLO 
contribution associated with the gg initial channel would ameliorate the situation and 
will be necessary to describe data at future high energy hadron colliders. 

Predictions are also presented in the paper for pN — > 77X at the center-of- 
mass energy 31.5 GeV, appropriate for the E706 fixed-target experiment at Fermilab. 
The large Qt and small A0 behavior of the kinematic distributions is dominated by 
the resummation of the gg initial state. Nonperturbative physics controls the Qt 
distribution, and neither the perturbative Sudakov nor the regular NLO contribution 
plays an important role, except in the very large Qt region where the event rate 
is small. For the E706 kinematics, the requirement Qt < Q is generally satisfied. 
Therefore, the approximate gg calculation presented in this work should be reliable. 

In this calculation, the incident partons are assumed to be collinear with the in- 
cident hadrons. A recurring question in the literature is the extent to which finite 



"intrinsic" &t may be required for a quantitative description of data [27, 17]. An 
important related issue is the proper theoretical specification of the intrinsic com- 



ponent [p8| . In the CSS resummation formalism, this physics is included by prop- 
erly parametrizing the nonperturbative function W NP {b), which can be measured in 
Drell-Yan, W, and Z production. Because photons participate directly in the hard 
scattering, because their momenta can be measured with greater precision than that 
of hadronic jets or heavy quarks, and because the 77 final state is a color singlet, the 
reaction pp — > 77X may serve as a particularly attractive laboratory for the under- 
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standing of the role of intrinsic transverse momentum. The agreement with data on 
the Qt distributions in Figs. £| and |7| is suggestive that the CSS formalism is adequate. 
However, the separate roles of gluon resummation and the assumed nonperturbative 
function in the successful description of the Qt distributions are not disentangled. 
In the non-perturbative function of Eq. dl2|) , the dependence on b (and, thus, the 
behavior of da/dQx at small Qt) is predicted to change with both Q and the values 
of the parton momentum fractions Xj. At fixed Q, dependence on the values of the 
X{ translates into dependence on the overall center-of-mass energy of the reaction. 
As data with greater statistics become available, it should be possible to verify these 
expectations. In combination with similar studies with data on massive lepton-pair 
production (the Drell-Yan process), it will be possible to determine whether the same 
non-perturbative function is applicable in the two cases, as is assumed in this paper. 

The diphoton data may allow a study of the nonperturbative as well as the per- 
turbative physics associated with multiple gluon radiation from the g/ttcm-initiated 
hard processes, which cannot be accessed from Drell-Yan, W ± , and Z data. With 
this knowledge, it may be possible to improve calculations of single photon produc- 
tion and other reactions sensitive to gluon-initiated subprocesses. In the D0 data 
analysis ||16|| , an asymmetric cut is applied on the transverse momenta (p T ) of the two 



photons in the diphoton event. This cut reduces the effect of multiple gluon radiation 
in the event. To make the best use of the data for probing the interesting multiple 
gluon dynamics predicted by the QCD theory, a symmetric p T cut should be applied. 
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Figure 1: Feynman diagrams representing the leading order and next-to-leading 
order contributions to photon pair production in hadron collisions. The shaded circles 
signify the production of long-distance fragmentation photons, which are described 
by the fragmentation function D^ q . 
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Figure 2: Comparison of the parton-level and Monte Carlo fragmentation contri- 
butions at the Tevatron. The upper and lower curves of the same type show the 
contribution before and after an isolation cut. The left figure shows the transverse 
momentum of the photon pair Qt- The right figure shows the light-cone momentum 
fraction carried by the fragmentation photon. 
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yy Resummation at Vs=1.8 TeV 




Figure 3: The predicted distribution for the invariant mass of the photon pair M 77 
from the resummed calculation compared to the CDF data, with the CDF cuts im- 
posed in the calculation. 
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yy Resummation at Vs=1.8 TeV 




Figure 4: The predicted distribution for the transverse momentum of the photon pair 
Qt from the resummed calculation compared to the CDF data, with the CDF cuts 
imposed in the calculation. 



33 




Figure 5: The predicted distribution for the difference between the azimuthal angles 
of the photons A(/> 77 from the resummed calculation compared to the CDF data, with 
the CDF cuts imposed in the calculation. 
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yy Resummation at Vs=1.8 TeV 




Figure 6: The predicted distribution for the invariant mass of the photon pair M 77 
from the resummed calculation compared to the D0 data, with the D0 cuts imposed 
in the calculation. 
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yy Resummation at Vs=1.8 TeV 




Figure 7: The predicted distribution for the transverse momentum of the photon 
pair Qt from the resummed calculation compared to the D0 data, with the D0 cuts 
imposed in the calculation. 
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yy Resummation at Vs=1.8 TeV 
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Figure 8: The predicted distribution for the difference between the azimuthal angles 
of the photons A0 77 from the resummed calculation compared to the D0 data, with 
the D0 cuts imposed in the calculation. 
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yy Resummation at Vs=1.8 TeV 




Figure 9: The predicted distribution for the invariant mass of the photon pair M 77 
from the resummed calculation. The additional cut Qt < M 77 has been applied to 
reduce the theoretical uncertainty. 
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yy Resummation at Vs=1.8 TeV 




Figure 10: The predicted distribution for the transverse momentum of the photon pair 
Qt from the resummed calculation. The additional cut Qt < M 77 has been applied 
to reduce the theoretical uncertainty. The lower solid curve shows the prediction of 
the pure NLO (fixed-order) calculation for the qq and qg subprocesses, but without 
fragmentation contributions. 
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yy Resummation at Vs=1.8 TeV 
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Figure 11: The predicted distribution for the difference between the azimuthal angles 
of the photons A0 77 from the resummed calculation. The additional cut Qt < M 77 
has been applied to reduce the theoretical uncertainty. 



40 



yy Resummation at Vs=31.5 GeV 
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yy Resummation at Vs=31.5 GeV 
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Figure 13: The predicted distribution for the transverse momentum of the photon 
pair Q T from the resummed calculation appropriate for pN — > 77X at v / S'=31.5 
GeV. 
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yy Resummation at Vs=31.5 GeV 
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Figure 14: The predicted distribution for the difference between the azimuthal angles 
of the photons A0 77 from the resummed calculation appropriate for pN — > 77X at 
^=31.5 GeV. 
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